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Ught with Fermi-IjAT and TeV observations of blazars 
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ABSTRACT 

We propose a method for setting upper limits to the extragalactic background 
light (EBL). Our method uses simultaneous Fermi-LAT and ground-based TeV 
observations of blazars and is based on the assumption that the intrinsic spectral 
energy distribution (SED) of TeV blazars lies below the extrapolation of the 
Fermi-LAT SED from GeV to TeV energies. By extrapolating the Fermi-LAT 
spectrum, which for TeV blazars is practically unattenuated by photon-photon 
pair production with EBL photons, a firm upper limit on the intrinsic SED at 
TeV energies is provided. The ratio of the extrapolated spectrum to the observed 
TeV spectrum provides upper limits to the optical depth for the propagation of 
the TeV photons due to pair production on the EBL, which in turn sets firm 
upper limits to EBL models. We demonstrate our method using simultaneous 
observations from Fermi-LAT and ground-based TeV telescopes of the blazars 
PKS 2155-304 and lES 1218+304, and show that high EBL density models are 
disfavored. We also discuss how our method can be optimized and how Fermi and 
X-ray monitoring observations of TeV blazars can guide future TeV campaigns, 
leading to potentially much stronger constraints on EBL models. 

Subject headings: diffuse radiation — galaxies: active — quasars: general — 
gamma rays: galaxies — quasars: individual (PKS 2155-304, lES 1218+304) 
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Introduction 



The EBL reflects the cosmologically important time-integrated history of hght produc- 
tion and re-processing in the Universe. For this reason, measuring its intensity is highly 
desirable. The two components of the EBL are dust emission peaking at A ~ lOO/xm and 
starlight peaking at ~ l/im. The actual level of the EBL is very difficult to measure, due to 
the dominance of foreground emission, mostly from iri terplanetary dust in our solar system 



(for reviews see iHauser fc Dwekll200ll : lKashlinskyll2005l ) , and the EBL level remains unknown 



within a factor of ~ few. Model-independent lower limits to the EBL, based on galaxy counts 
(e.g. Dole et al. 2006, Bethermin et al. 2010), are the only strict lower limits on the EBL 
to date. Modeling the EBL can lead to definite prediction, however uncertainties in the star 
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proposed a new method based on detecting as GeV emission EBL radiation that has been 
inverse Compton-scattered by relativistic electrons in the lobes of nearby radio galaxies such 
as Fornax A. 

The EBL in the 1-10 /xm range can in principle be obtained by using the TeV blazars 
as background light sources and modeling its attenuation due to pair production with EBL 
photons: by assuming that the intrinsic TeV spectrum is known from modeling of the broad- 
band blazar SED, one can deri ye the mid-IR EBL by comparing the observed to the pre- 



sumed intrinsic spectrum ( e.g., IStecker. de Jager. fc Salamonlll992l : IStanev fc Franceschini 



19981 ; [Renault et al. 1 1200 ll ) . Because it is not possible to determine with confidence the 
intrinsic TeV spectrum, a variation of this method has been proposed that sets limits 
on the EBL by assuming that the intrinsic blazar TeV spectrum cannot be arbitrarily 
hard: from simple shock acceleration argum ents one would not ex pect an intrinsic TeV 



photon index harder than TrpeV = 1-5 (e.g., lAharonian et al.l 120061 ). Detailed shock ac 



celeration simulations, however, indi cate that harder VHE spectral indices may be possible 
( Stecker. Baring. &: Summerlinll2007l ). A large lower electron Lorentz factor (IKatarzyhski et al 
20061 ). and Compton scattering of the cosmic rnicrow ave background in an extended jet may 
also lead to hard TeV spectra (iBottcher et al.l 12008). These considerations significantly re- 
lax th e EBL limits derived by assuming T-TeV > 1-5 (iMazin fc Rauell2007l : lFinke fc Razzaque 
2009h . 



Methods that constrain the EBL through purely spectral arguments are free of the 
uncertainties of adopting a pa rticular physical mo d el. Su ch methods, based solely on TeV 
data, have been proposed by iDwek fc Krennrich I (120051 ). who considered unnatural TeV 
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SEDs that exhibit an exponential increase at their high energy end and by lSchrodter I (120051 ) 
who assumed that all TeV blazars in flaring states have TeV spectra with the same maximum 
intrinsic photon index Txev = 1-8. Here we present a spectral method for obtaining upper 
limits to the EBL energy density that makes use of simultaneous LAT and TeV observations. 
In §2] we describe our method and apply it to recent simultaneous LAT and TeV observations 
observations of PKS 2155-304, in ^we discuss how we can produce stronger constraints on 
the EBL and demonstrate this using simultaneous GeV/TeV observations of lES 1218+304, 
and in SHlwe conclude. 



2. An upper limit for the intrinsic TeV SED of TeV blazars 



Our current observational understanding of blazars points toward a SED that consists 
of two spectral components or more colloquially "bumps" . The first bump peaks at IR to X- 
ray energies and it is almost certainly synchrotron emission from a population of relativistic 
electrons in a partially ordered magnetic field. The second bump peaks at MeV to multi- 
GeV energies and is thought to be due to inyerse C ompton scatter ing of synchrotron photo ns 



(IBand &: Grindlaylll986l: iBloom fc Marscherlll996[) . a dust y torus (IBlazeiowski et al.ll20001). a 



broad-line region ( ISikora et al.lll994j ). or an accretion disk ( iDermer. Schlickeiser. &: Mastichiadis 
19921 ). For the blazars that have been discovered so far, this second bump ranges in power 
between 0.1 — 100 times the power of the synchrotron bump. It is also possible that the 
high energy bump is a result of emis sion from hadrons co-acce lerated with the jet, which 
can radiate by hadronic synch rotron (iMiicke fc Protherodl200ll ) or photomeson production 
(IMannheim fc Biermanrull992l ) and the resulting cascades. Protons may a lso convert to neu- 
trons , escape the relativistic jet, then convert back to protons and radiate (lAtoyan fc Dermer 
20031 ). So far, no observational evidence for a third, higher energy bump has been found and 
we consider it unlikely that a third SED component at TeV energies will be more powerful 
than the extrapolation of the GeV SED at TeV energies. 

Our working assumption is motivated by these considerations and is weaker than the 
statement that there is no third SED bump: we assume that the extrapolation of the Fermi- 
detected SED at TeV energies is higher than the intrinsic TeV SED of the source. This 
means that the ratio of the actually observed TeV fiux f^^obs to the extrapolated one fe,ext 
at any given TeV energy e, provides a strict upper limit r^^max to the EBL-induced pair 
absorption optical depth at this energy. 



MUext/ fe ,obs I 



This optical depth can then be compared to the optical depth calculated for the various EBL 
models. The models for which the optical depth is greater that ^ax are then excluded. 
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An important consideration for our method stems from the mismatch of the relatively 
long time required to obtain a sufficiently good SED of typical TeV blazars with Fermi (no 
less than ~ a day) to the significantly smal ler variability time ob served by TeV telescopes 
(in some cases lasting less than an hour, e.g. lAharonian et al.ll2007l ). It is clear that, because 
we want to compare simultaneous spectra, the TeV fluxes observed during the rapid TeV 
variability events seen in TeV sources should not be used together with the Fermi SED 
derived after integrating for days. 



2.1. Application to PKS 2155-304 



The blazar PKS 2155-304, a high peak frequency BL Lac object at redshift z = 0.116, 
was the target of a multiwavelength campaign in late August and early September 2008, 
which included observations by the Fermi Gamrna- Ray Space Telescope and the HESS at- 
mospheric Cherenkov telescope (lAharonian et al.l l2009bl ) . It is these observations to which 
we now turn ou r attent ion. The 7-ray SED from this campaign can be seen in Fig. [T] 
Aharonian et al.l (l2009b( ) found that the Fermi SED derived from data between MJD 54704- 
54715, the period of the TeV observations, can be described by a simple power law with 
spectral index F = 1.81 ± 0.14. A longer train of Fermi observations between MJD 54682- 
54743 exhibits a similar GeV state. If one includes these data, the exposure increases by 
a factor of ~ 3.6 and the preferred GeV spectr um is now a broken pow er law, whose high 
energy part has a photon index F = 1.96 ±0.11 ( lAharonian et al.ll2009bl ). While a conserva- 
tive estimate of the photon index used for obtaining an upper limit to the TeV flux is that 
derived from the simultaneous observations only, the lack of strong variability in the GeV 
(as well as the TeV) regime, suggests that it is reasonable to adopt the high energy photon 
index better description. 

The single power-law, with spectral index F = 1.81 ± 0.14 and high energy power-law 
flt, F = 1.96 ± 0.11 of the broken power-law are shown as "bow tie" error plots in Fig. [H 
extrapolated to the HESS energy range. These extrapolations are used as upper limits to the 
intrinsic flux of the TeV SED, unabsorbed by the EBL. The upper limit on r-y-y is calculated 
from equation ([1]) for these two extrapolations, and the results are shown in Fig. |2] along with 
the absorption optical depth predictions of several EBL models. As expected, the steeper 
GeV index provides stronger constraints to the EBL models. Note that the fast evolution 
model of Stecker et al. 2006 lies below the la limit on r^^ at the highest TeV energy for 
F = 1.96 ± 0.11 and, in fact, is inconsistent with it at the lAa level. 
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3. How to produce stronger constraints on the EBL 

We now discuss liow tlie yield of our method can be maximized. We base our discussion 
on the fact that Fermi is continuously monitoring the entire sky, including the extragalactic 
TeV sources. To facilitate our discussion consider a source that has a LAT spectral index 
aLAT and a TeV spectral index axev (spectral indexes a are connected to photon indexes T 
through the relation a = T — 1). Let us further assume that the transition between the two 
power laws takes place at an energy ei ~ 100 GeV, practically the energy border between 
LAT and ground based TeV telescopes. Then at an energy £2 within the TeV regime, the 
flux one would expect by extrapolating the LAT spectrum is fe2,ext = /ei(e2/ei)~°^^^, while 
the flux fe2,obs "we actually observe at energy €2 is fe2,obs = /ei(e2/ei)~"^''^, where /^j is the 
flux at energy ei. This means that the maximum pair production optical depth r^j.max that 
we can infer from observations of this source is 

Te2,max = {ctTeV " CXlAt) In {^~^ ' (2) 

Note that since T^2,max increases with e2, the higher the TeV energy in which a source is 
detected, the stronger the constraints will be. As expected, the smaller the spectral break 
Aa = aTeV — OiLAT of a source is, the stronger the resulting constraint on T^2,max will be. 
Therefore, at a given redshift z, the most promising sources are those that exhibit the smallest 
Aa. If a source varies, states with the smallest Aa will provide the strongest constraints. 
We discuss this further below. 



3.1. The benefits of high TeV states 



A general pattern observed in both the low and high energy bumps of blazar SEDs 
is that variability events us ually manifest the r nselves with a hardenin g of the high energy 
part/tail of the bump (e.g., lAlbert et al.l 120071 : lAharonian et al.ll2009al ). with the amplitude 
of variability decreasing at progressively lower energies. It is, therefore, natural to anticipate 
that when a TeV blazar is in a high state, the spectrum will harden and rise in amplitude 
more at TeV than at GeV energies (variability can explain the fact that in some variable 
sources like MRK 421 or S5 716+714 the extrap olation of the GeV SED of a given non- 
flaring epoch to TeV energies (lAbdo et al.ll2009bl ) can be below the actually detected TeV 



emission at flaring states). This means that in a flaring state the method will provide more 
severe limits on EBL models. 



As a potential example of the benefits of catching a source at a high state, we return 
to the case of PKS 2155-304, whose TeV flux in the August- September 2008 Fermi-HESS 
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campaign was close to the lowest archival TeV data. The source exhibited a spectral break 



Ag = 1.38 {a^AT = 0.96, axev = 2.34). In a pre- Fermi 2006 campaign (lAharonian et al. 



2009al ). the TeV spectral index during flaring states hardened up to a^ev = 1-6. If a 
similarly hard TeV spectrum is recorded simultaneously with Fermi observations and if the 
Fermi spectrum does not harden significantly, then Aa = 0.74, which would reduce Tmax by 
a factor of ~ 2, posing important constraints on the EBL (all the r^^ upper limits in Figure 
[2] would shift downward by a factor of ~ 2). 

The ideal high states should last long enough to provide us with a solid determination 
of the Fermi fiux and photon index. In addition, the source should remain relatively steady 
while in the high fiux state, so that the average Fermi and TeV measurements are a good 
representation of the source during this high but not rapidly variable state. A high level of 
GeV fiux is also very useful, because if there is curvature in the Fermi band, it is possible 
that a broken power law fit will be better than a simple power law and one can then use the 
higher energy steeper part of the Fermi SED to extrapolate to TeV energies and produce 
lower values of Tmax (as we did for PKS 2155-304). 



3.2. Sources at diflTerent redshifts, the case of lES 1218-^304 

Even if there are sources described by pure power laws in the entire GeV - TeV regime 
without the need for a break, their observed spectra would be imprinted with a break solely 
attributed to EBL absorption. Such breaks would increase with increasing redshift. In the 
more pragmatic case that there is a distribution of intrinsic breaks, then, as the blazar 
sequence (Fossati et al. 1998) suggests, the breaks should become stronger for the more 
powerful sources seen at higher redshifts, because the peak frequency of the high energy 
component shifts to lower energies as the source power increases. The increase of Aa with 
redshift would thus be a convolution of two effects: the blazar-sequence-like shift of the 
TeV spectrum to steeper values (which is intrinsic) and the increase of EBL absorption with 
distance. 

If Aa is dominated by the former, our method would derive its strongest EBL constraints 
from nearby sources. If, however, the intrinsic increase of the break with redshift is small or 
negligible, sources at higher redsh ifts, even with steep TeV spectra (the most distant example 



is 3C 279 discovered by MAGIC flAlbert et all 120081 ) at z = 0.538 with TtcV = 4.11 ± 0.68) 



would provide us with strong constraints on the EBL. Because we do not actually know 
the intrinsic breaks of the TeV sources, the question that naturally arises is under what 
conditions higher redshift sources with steep TeV spectra can provide as useful constraints 
on the EBL as their nearby siblings that have significantly harder TeV spectra. 



- 7- 



The important thing to notice here is that the relevant quantity that constrains the 
EBL is not r^ax, but T^ax/dL, where di is the luminosity distance of the source: sources 
with the same value of T^ax/dL will provide the same constraints on the EBL (assuming 
that the evolution in the EBL energy density itself with redshift is negligible for small 
redshifts). This means that, for fixed values of ei and €2, a particular constraint on the 
EBL is satisfied for sources with a fixed value of Aa/di. More distant sources with larger 
breaks Aa will provide the same constraints as more nearby sources, as long as they have 
the same Aa/d^. As an example, let us use as a reference the state of PKS 2155-304 we 
discussed above, with Aa = 1.38 and cIl = 533.4 Mpc. A source at the distance of Mrk 
421 {z = 0.034, = 143 Mpc) will provide the same constraints on the EBL if it exhibits 
Aa = 1.38 X 143/533.4 = 0.37. Any break greater than Aa = 0.37 will not be as constraining 
as the stronger break of the more distant PKS 2155-304. 



^Going to higher redshifts, the more distant TeV BL Lac lES 1218+304 (lAlbert et al. 

20061 ) {z = 0.182, di = 873.5 Mpc) would provide the same constraints on the EBL if it 
exhibits Aa = 1.38 x 873.5/533.4 = 2.26. Breaks gentler than that will produce stronger 
constraints on the EBL than those derived from the state of PKS 2155-304 we studied here. 
Recently, it has been reported that PKS 1218+304 has been observed by VERITAS from 
2008 December 29 to 2009 April 23 (Acciari et al. 2010). This is quasi-simultaneous to the 
LAT observation of the source in the 11 month Fermi-LAT catalog (A. Abdo et al. 2010 in 
preperation)0, which continuously exposed PKS 1218+304 from 2008 August 4 to 2009 July 
4. These LAT and VERITAS spectra give Aa = 1.37 indicating they can provide a stronger 
constraint on EBL models. 

Following the same procedure as in the case of PKS 2155-304, we plot in Figure [3] the 
7-ray SED, together with the 'bow tie' extrapolation of the GeV spectrum to VERITAS 
energies. As before, the extrapolation is used as an upper limit to the intrinsic flux of the 
TeV SED and the upper limit on r^^ is calculated from equation ([T]) and plotted in Fig. H] 
along with the absorption o ptical depth pr e dictio ns of several EBL models. As can be seen. 



the fast evolution model of lStecker et all n200a) lies below the la upper limit on r-y-y for 



all but the two lowest TeV energies, and below the 3a level at the highest TeV energy. In 
particular, the hi ghest point 1.8 T eV) is inconsistent with Stecker's fast e v olutio n model 
at 4.7(7. Also the lStecker et al\ feOOq ) baseline model, and the lKneiske et all l(2004) best fit 



model are inconsistent at the 2.6a level and 2.9a level with Tmax at the highest TeV energy. 



^The catalog is published online at http : //f ermi . gsf c ■nasa.gov/ssc/data/access/lat/lyr_catalog/| 
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Conclusions 



We presented a simple, model-independent method for setting upper limits to the EBL. 
Our method is based on the assumption that the level of the intrinsic TeV emission of blazars 
is below the extrapolation of the LAT SED to TeV energies. We applied our method to PKS 
2155-304 and lES 1218+304, the only two TeV blazars of known redshift for which published 
simultaneous LAT-TEV observations currently exist. Even with these first applications of 
our method, the highest level EBL models are disfavored. Future LAT- TeV simultaneous 
observations hold the promise of pushing the EBL upper limits much further. We argued 
that it is important to devote TeV time not only to nearby sources, but also to more distant 
sources, hoping in both cases to observe high TeV states that for a given source will exhibit 
smaller GeV to TeV spectral break Aa, and will, therefore, produce stronger constraints. 

Because the most difficult observations to obtain are in the TeV band, requests to 
monitor particular sources with TeV facilities can be triggered from Fermi observations of 
high states of TeV sources. A complementary approach is X-ray monitoring of the TeV 
blazars. In this case, because the X-ray emission of most TeV blazars is the high energy tail 
of the sync hrotron component, high X-ray states are, in general, a good proxy for high TeV 
states (e.g. lAharonian et al.l (l2009aj)). Such X-ray monitoring holds the promise of catching 
states in which, while the GeV emission does not increase substantially, the TeV emission 
does. 

We note here that our method assumes that the entire spectral break from the LAT to 
the TeV bands is due to EBL pair production absorption. This is an extreme assumption 
and it is highly probable that a substantial fraction of the break is intrinsic to the source. 
This in turn means that the actual level of the EBL may be significantly lower than the 
upper hmits produced by our method. It would be very exciting and possibly hinting to new 
physics (e.g., lAmelino-Camelia fc Piran 1 120011 ) if the lowest collective values of Tmax that 
our method will produce, cha l lenge the lower leyel on the EBL inferred by galaxy counts 



(e.g., iMadau fc Pozzettil I2OOOI : iFazio et al. 1 12004 iBethermin et al. II2OIOI ). We anticipate 
that current and upcoming TeV- GeV blazar monitoring campaigns will provide plenty of 
opportunity for applying our method. 
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Fig. 1. — The Fermi data (black diamonds) and the extrapolated Fermi SED for the two 
cases described in the text, together with the HESS data, (red diamonds) for PKS 2155-304. 
The solid lines depict the energy range actually covered by LAT, while the broken lines 
depict the extrapolation of the LAT SED in the TeV regime. 
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Fig. 2. — The pair production optical depth along the line of sight to PKS 2155-304 for a 
range of EBL models, and the constraints our method imposes on the actual optical depth for 
the two cases for the GeV slope described in the text. Empty and solid triangles correspond 
to la and 3a upper limits on r^^. 
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Fig. 3. — The Fermi-LAT SED {Fermi 1 year on-line LAT catalog; solid lines) extrapolated 
to TeV energies (broken lines), together with the VERITAS data (red diamonds) for lES 
1218+304. 
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Fig. 4.— Same as in Fig. El but for lES 1218+304. 



